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The conformal anomalies and functional determinants of the Branson- 
GJMS operators, P 2k , on the (i-dimensional sphere are evaluated in 
explicit terms for any d and k such that k < d/2 (if d is even). The 
determinants are given in terms of multiple gamma functions and a 
rational multiplicative anomaly, which vanishes for odd d. Taking the 
mode system on the sphere as the union of Neumann and Dirichlet 
ones on the hemisphere is a basic part of the method and leads to a 
heuristic explanation of the non-existence of 'super-critical' operators, 
2k > d for even d. Significant use is made of the Barnes zeta function. 
The results are given in terms of ratios of determinants of operators on 
a (d + 1) -dimensional bulk dual sphere. For odd dimensions, the log 
determinant is written in terms of multiple sine functions and agree- 
ment is found with holographic computations, yielding an integral over 
a Plancherel measure. The N-D determinant ratio is also found ex- 
plicitly for even dimensions. Ehrhart polynomials are encountered. 
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1. Introduction. 

The conformally covariant higher derivative GJMS operators, P 2 k, [1], have 
been the subject of a certain amount of activity in the general area of conformal 
geometry e.g. Juhl, [2], Gover, [3]. The paper of Diaz, [4], contains some physical 
applications and a useful survey. 

Branson has computed the ratio of the (log) determinants of such operators in 
conformally related spaces, and has defined the notion of Q-curvature to be of use 
in this endeavour. The method is the standard one, well known to physicists since 
the late 1970s, of integrating the conformal anomaly, which is the constant term 
in the heat-kernel expansion, or the value of the ^-function at zero, and is a local 
quantity easily computed. 

If one wishes to know the actual value of the determinant using this Jacobian, 
or cocycle function, it is necessary to evaluate it in some fiducial metric where 
the eigenvalue problem is sufficiently tractable, in one way or another. Tradition- 
ally, spheres and generalised cylinders have been used for this purpose, in other 
situations. This is a more difficult problem as it leads to non-local quantities, in 
particular the derivative of the ^-function at zero, in ^-function regularisation. 

Branson, [5], constructed the GJMS conformally covariant higher derivative 
operators in the special case of the (round) sphere as an explicit product, and I 
intend here to use direct, simple-minded spectral techniques to work out the value 
of the ^-function of P 2 fc a ^ argument (the conformal anomaly) and of its derivative 
at 0, which, as mentioned, determines the determinant. 

For comparison, Guillarmou, [6], has already computed the determinant of _P 2 fc? 
holographically, on the conformal infinity of a factored hyperbolic space. Trivial 
factoring gives the case of a bounding sphere. See also Diaz, [4]. Diaz and Dorn, [7], 
gave an independent conformal field theory computation of an equivalent quantity 
on this boundary using dimensional regularisation, again in a holographic setting. 



2. The Branson GJMS operators 



I denote by Y d the (non-negative) conformally invariant Laplacian, 
often referred to, in mathematical works, as the Yamabe operator. 
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Branson's construction of P 2k , [5], in the special case of the (round) <i-sphere 
is a simple product 2 , k G Z + , 

fc-i 

^*=n( y d-^' +1 )) 

(i) 

fc-i v > 

= J] + - j - 1/2) (0^ + 1/4 + j + 1/2) . 

There are not many derivations of this factorisation. Branson uses Lie group 
theory to construct an object which has the required properties and then invokes 
a uniqueness argument. Amongst the few other derivations, Graham's, [8], is the 
most direct, using stereographic projection from flat space. Morpurgo gives a similar 
discussion earlier, [9], and makes contact with Branson's method. 

I assume, in a physicist's non-rigorous way, that the eigenvalue problem for 
P 2k is solved by that for Y d , and this is classically known. It has been shown, 
for example, with a certain amount of effort, that the P 2k are self-adjoint. Some 
function analytic discussion can be found in [9]. 

The reason for writing P 2k in the above way is because the eigenvalues of 
Y d + 1/4 are perfect squares. Branson, [5], also gives this form. His notation is 
B = y/Y d + 1/4. 

Furthermore, I am going to obtain the sphere expressions by adding the hemi- 
spherical Dirichlet and Neumann ones, as in previous work, [10-12]. The eigenvalues 
of y/Y d + 1/4 then take the general linear form m . u> + a where m and u) are d— 
dimensional vectors with m i G Z + and (here) uj i = 1, Vz. 3 This is the way the 
eigenvalues immediately appear in the classic, and ancient, separation of variables, 
the degeneracies arising from coincidences. It is better not to do the combinatorics 
which give them explicitly in terms of factorials. 

The constant, a, takes the values a = a N = (d — l)/2 and a = a D = (d + l)/2 
for, respectively, Neumann (N) and Dirichlet (D) conditions on the rim of the 
hemisphere. I write out the eigenvalues of P 2k , 

k-i 

A fc (m) = ] \ (m . uj + a - j - 1/2) (m . w + a + j + 1/2)) , (2) 

3=0 

2 I make a different sign choice to others. 

3 I retain u> for a while because the expressions have validity when the sphere is replaced by a 
quotient. 
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and remark that the first factor has a zero in the Neumann case at the upper 
limit first, as k increases, when k = d/2. (I am assuming initially that d is even). 
Increasing k to d/2+1 and over, gives zeros for N, but also, more vitally, for Dirichlet 
conditions, which is impossible 4 , and so we must restrict k to k < d/2. The available 
general proofs of this depend upon the existence of appropriate tensors. 

In the calculations, the critical case, k = d/2, must be treated separately. It 
corresponds, roughly, to the minimally coupled scalar. I now turn to the construc- 
tion of some spectral invariants. 



3. The zeta function and its value at zero 

In the usual way, I define the relevant (^-function as the continuation of the 



sum, 



oo ^ 

Z d {a,a,k)= E^TyTJ 

oo fc-1 „ (3) 



En 



m=U j 



' A = A ((m. W + o)2-0- + l/2) 2 )- 



and am required to find Z d (0,a, k) and Z' d (0,a, k). 5 The evaluation at s = is 
the condition that allows an explicit solution along the lines followed in [11]. This 
reference deals with the lowest case, k = 1. In order to make comparison simpler I 
rewrite (3) 



Z d (s, a, k) = 77 7" v2 Ti —\~2 i - ^' ( 4 ) 

^ ((m . w + a) 2 - a 2 ) . . . ((m . w + a) 2 - a^J 

where = j • + 1/2. In [11] it was shown that, in order to deduce Z d (0, a, k) and 
Z' d (0, a, k), it was sufficient to expand in powers of a and I will employ the same 
tactic here, for every factor. The procedure is adequately illustrated by the k = 2, 



4 The argument is an extension of the one familiar in potential theory and depends on the 
alternating signs of the powers of the Laplacian in P 2 fe to gi ye a sum °f squares on application 
of Green's theorem. If 2k > d, we are led to operators that have a Dirichlet mode on the 
hemisphere. This is impossible and so, while one can write the operators down, they have to be 
multiplied by zero, which is always an option. 

5 Another task would be to determine the singularity structure. 
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Paneitz operator case, 



Z d (s,a,2) = 



p ((m. uj + a) 2 — a 2 ) ((m . u> + a) 2 — a' 2 )' 

E 



m- 

v 2r , 2r ,s( S +l)...(s + r-l) s(s+l)...(s + r'-l) (5) 



a a 



r,r'=0 

x C d (4s + 2r + 2r / ,a|w), 

with a = 1/2, a' = 3/2 and where the mode sum has now been performed to give a 
Barnes ^-function as in [11], and elsewhere. This expression does not constitute a 
complete continuation. It does not give the poles, for example, but the values when 
s is a non-positive integer can be obtained fairly easily. 
The general definition of the Barnes ^-function is, 



iT(l — s) f exp(— ar)(— r 



,s-l 



q g ,a|uQ= Zr(1 ~ g) ( dr- 

^ Jl nti (1 - exp(- Wi r)) 

1 ( 6 ) 
= > r- , Re s > d , 

m=0 v ' 

where I refer to the components, u> i , of the rf-vector, u>, as the degrees or parameters. 
For simplicity, I assume that the uo i are real and positive. For the hemisphere they 
are all equal to 1, as given above. If a is zero, the origin, m = 0, is to be excluded. 
The contour, L, is the standard Hankel one. 

If s is set to zero in (3) one has to use the fact that the Barnes (^-function has 
poles, 

Q(s + r, a [ u?) -» ^ + Rr(d, a) as s ->■ 0, (7) 
where 1 < r < d. Then very easily, 

Z d (0, a, 2) = C d (0, a\u>) + -J2 " ^(4 a) , 

r=l 



where u = d/2 if <i is even and u = (d — l)/2 if <i is odd. 

The residues are known in terms of generalised Bernoulli polynomials but they 
are not needed because it was shown in [11] that there is an identity, 



r=l 



" a 2r 

CM a + « |w) + Cd(0, a - a |w) - 2C d (0, a |w) = J] — N 2r (d, a) , (8) 
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and then, 

Z d (0, a, 2) = 1 (C d (0, a + a |w) + C d (0, a - a |w) + &(0, a + a' |w) + C d (0, a - a' |w)) , 

which is the average of the regularised dimensions of the linear operators occurring 
in (2). 

It is clear that the same structure will hold for all levels, and so, for the sub- 
critical cases, 

^ fc-i 

Z d (0, a,k) = —J2 (C d (0, a + aj ) + C d (0, a - a J) , k < d/2 . (9) 

3=0 

To obtain the hemisphere expression, the appropriate a coefficient has to be 
used. For the Dirichlet case, (9) holds for the critical k = d/2 case also, while for 
Neumann conditions, the existence of a zero mode means that 1 has to be subtracted 
from the right-hand side. 

Barnes, [13] gives the required values of the ^-function, in particular, 

Cd(0,a|«) = ^flJV|«) 
so that from (9), if d is even, 

fc— l 

Z M a D , k) = ^J2 ( B d d) ( d / 2 +J + 1) + ~ J)) 



fc-i 



= ^ E [Bf{d/2-j - 1) + Bf\d/2 + j)) 
= Z d (0, a N , k) , 



(10) 



where now I have dropped the label u> and used a symmetry of the Bernoulli poly- 
nomials, which shows that the N and D quantities are equal, for k < d/2. Adding 
these two then just doubles the expression and so, on the full i.e. periodic sphere, 



i k ~ 1 

Z d (0,k) 



= TdiH ( B( d d \d/2+j + 1) + Bf\d/2 + j)) - 5 kM , (11) 

sphere ft M - J=0 



which is the expression quoted in an earlier work, [14], where some further comments 
can be found. 
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The expression (11) can be re-expressed in various ways. Use of the integral 
relation 

B { ^- l \x) = [ dtBi n \x + t), 
Jo 

and the product form, (e.g. [15] p. 186, [16], §8), 

B { f +1 \x\ 1) = (x - l)(x - 2) . . . (x - d) , 
yields the compact form, 

2(-i) d / 2 rk d / 2 ~ 1 



Z d (0,k) 



+ 5k OA = 

sphere 



which is identical to a formula deduced by Diaz, [4], on a holographic argument. 

If d is odd, the (anti-) symmetry of the Bernoulli polynomials shows that 
Z d (0, a D , k) = —Z d (0, a N , k) and so the conformal anomaly on odd spheres is zero, 
as it must be. 

4. The determinants in subcritical cases 

It is necessary to find the derivative Z' d (0,a,k). I will again use the k = 2 
expression (5) to exemplify the calculation. 

The sums over r and r' are split up to give 

s(s + 1) . . . (s + r - 1) 



Z d (s, a, 2) =C d (4s, a) + X> 2r + a^) ^ ^ - ^ C d (As + 2r, a) 

, 2r /s(s + l)...(s + r-l) s(s + l)...(s + r'-l) (12) 



oo 



+ £ 



2r . 

a a 



r,r' = l 

x Q(4s + 2r + 2r',a) . 

The calculation of the derivative of the first two terms proceeds as in [11]. The 
difference is the sum (a 2r + a' 2r ) in place of a 2r and the argument of 4s instead of 
2s. Carrying through the algebra is straightforward and yields for the derivative at 
zero 

Z' d (0, a, 2) = Cd(0, a + a) + &(0, a - a) + &(0, a + a') + &(0, a " «') 

-it ^^ H 2 (r)N 2r (d,a)+X : (13) 



where H k is defined by, 

H k (r)=H(2r- 



in terms of the harmonic series 




X is the contribution of the third term in (12), now to be calculated. This 
term is the novelty. 

The pole structure, (7), and the existence of a factor of s 2 in the summand 
quickly shows that the remainder term, R, can be ignored and also that the sum 
over the region where the Barnes ^-function has no poles can be dropped. From the 
argument of the (^-function this region is r + r' > u. The sum is therefore restricted 
to a finite set of terms, which is a big calculational simplification. 

Taking the elementary derivative, and setting s to gives, combined with (13), 
the complete derivative at zero, 



The first four terms are to be expected from the complete linear factorisation of 
the eigenvalues, (2). The remaining ones are evidence of multiplicative anomalies 6 . 
The fifth term is the anomaly associated with the factorisation of each a 2 bracket, 
of which there are two. I term this a first order anomaly. The last term is a mixing, 
or cross term, corresponding to an interaction between the different a 2 factors. I 
call this a second order anomaly. 

If one raises the level to k by adding further factors to (5) one sees, by spitting 
the sums, that there are k contributions of first order and k(k — 1)/2 of second. Fur- 
ther, the increasing powers of s in the summand show that there are no third order 
contributions. The calculation of the determinants for the P 2k is then straightfor- 
ward and a matter of bookkeeping. 

6 In odd dimensions there is no anomaly or, rather, the N and D hemisphere anomalies cancel. 



Z' d (0, a, 2) = Cd(0, a + a) + Cd(0, a - a) + ^(0, a + a') + &(0, a - a') 




(14) 
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The derivatives of the Barnes function are formally expressed in terms of mul- 
tiple gamma functions, by definition of the latter, and I write the general log deter- 
minant (p d is a normalising modulus, [13], and is independent of a), 

Z' d (0, a, k) = log n F d( a + a j) F > - a i)) + M ^ a ' fc ) ' ( 15 ) 

where M(<i, a, fc) reflects the multiplicative anomaly and is composed of first and 
second order terms, M(d, a, k) = M x (d, a, k) + M 2 (d, a, k), given by, 

u / k—1 \ 

M x (d,a,k) = -J2[Yl a T )- H k( r ) N 2r (d, a) 

r=l ^ j=0 ' T 

u u-r i k-1 , . 

M 2 (d,a,k)=-J2 E E«N/ J 7^*-!)+-^--^ WV 2Pi+2r (d,a), 

(16) 

respectively. 



5. Computations 

Expression (15), with (16), gives the derivative of the (function for P 2k on 
the hemisphere for N and D conditions after substitution of the appropriate values 
for the a and the a. Adding these gives the full sphere value which is the one most 
frequently requested. 

The first part of (15) is difficult to find (numerically), being transcendental. 
The second part is easy, being rational. I only need to give the residues, 

= (r - l)\(d - r)\ B * }M ~ a) ' 
again in terms of generalised Bernoulli polynomials, [15,13]. These are easily cal- 
culated from recursion and the symmetry shown in (17) shows that the N and D 
residues are equal, N r (d,a N ) = N r {d,a D ). M is well adapted for evaluation by 
machine algebra. 

It might appear that the transcendental part of (15) becomes more involved for 
increasing k. I will now show this is not so, at least formally, by using the relation, 
([13]), 

T d (a) T d+1 (a) 



p d T d+1 (a + l) 



(18) 
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to effect cancellations which also eliminate the moduli p d . Adding the N and D 
contributions, and substituting the explicit a J7 the relation (18) allows staggered 
numerator-denominator cancellations which result in the compact formula, 



Z' d (0,k) 



^ T d+1 {d/2-k)T d+1 {d/2-k + l) 

sphere ^ d+1 {d/2 + k) T d+1 {d/2 + k + 1) 



dropping the a N argument by defining M(d, k) — M(d, a N , k). 

I do not intend to evaluate the first term numerically nor to take it much 
further formally for even d (see section 10). An expression for logT n as a sum of 
derivatives of the Hurwitz (-function can be found, e.g. by using the expansion in 
[16] §9, and elsewhere. 

However, I do present some numbers for the multiplicative anomaly part ob- 
tained from (16) in the form of polynomials in k, for a few d. There is no anomaly 
for k = 1/2, i.e. for just one linear factor, and so I list the polynomials M(d, k) = 
M(d, fc)/(l-2fc), 

M(4, k) = ^gpC 1 + 2£;)(132£; 3 - 36k 2 - 137k + 31) 

M(6, k) = (13152/c 6 + 2256/c 5 - 79344/c 4 - 22032/c 3 

v ' 29030400 v 

+ 76330/c 2 + 26195/c - 5805) , 

and give the values for the Paneitz operator for different dimensions, M(4, 2) = 
-659/1152 w -0.57205, M(6,2) = 5141/276480 w 0.01859 and M(8, 2) = 
-15666659/6502809600 w -0.00241. 



6. The determinant in the critical case 



When k = d/2 there is a zero mode in the Neumann problem and it is necessary 
to remove this from the ^-function, i.e. the entire m = term has to be omitted 
in (3), and hence in all the arising Barnes ^-functions. 

The first thing to say is that the multiplicative anomaly contribution is un- 
altered because removing one term does not affect the analytic structure of the 
^-functions and, although it does affect the value at zero, this cancels in the rele- 
vant identity, (8). 

Referring again to the k = 2 example, it is the last C'(0) term, with a N = a' , 
in (14) that gives the trouble, solved by replacing all the Barnes (-functions by 
subtracted ones, called ( d in [11]. This can be easily allowed for in the other ('(CI) 
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parts by adding and subtracting their m = terms, producing a correction. For 
the offending term, however, I have to use the limit, [13], 



iimCdtM = -loge-logp d . 

e— >0 



These adjustments mean that (19) is replaced by 



Z'Md/2) 



sphere 



= - log 



(d-l)\T d+1 (d)T d+1 (d+l) 

r d+1 (i) 

+ 2M(d, d/2) . 



(20) 



The cancellation pattern has altered, slightly, and the correction has been 
added. I also note that p d = r d+1 (l). 

7. Explicit evaluation in the odd case 

As has been mentioned, if d is odd, the anti-symmetry of the appropriate 
Bernoulli polynomials means that the multiplicative anomalies for N and D condi- 
tions have opposite signs, and so cancel for the complete sphere. Furthermore, P 2k 
has no zero modes 7 , in either case, and so the result for the log determinant is as 
in (19) except there is no final M term, and all k are allowed. The only difference 
is that the existence of negative eigenvalues for k > d/2 introduces an overall sign 
factor, (-l)fl( A <°) (see section 8). 

It is well known that the odd and even <i-spheres form two, somewhat distinct, 
families. In the former, for conformal coupling in d + 1 dimensions, the heat- 
kernel expansion terminates 8 . The heat-kernels can be related by a derivative 
which increases d by 2 and the simpler odd dimension properties can be traced to 
the facts that S 1 and S 3 are group manifolds. 

Another consequence is that C'(0) ( equivalent ly, the effective action) can be 
computed in terms of 'elementary' functions and I will perform the calculation that 
illustrates this using (19), which I repeat, 



Z' d (0,k) 



sphere 




(21) 



= log Sin d+1 (d/2 + k) - log Sin d+1 (d/2 - k) , d odd , 



There is no analogue of minimal coupling in the class of P 2 fe operators for odd d. 
In the dual hyperbolic space this expansion is exact. 
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where I have employed Kurokawa's multiple sine function S d (x) = Sm d (z) (see the 
appendix). For computation purposes, it is best to relate this to the the multiple 
cotangent function (see the appendix) and write (21) as an integral. This will also 
enable me to compare with other formulae. Thus, 



Z' d (0,k) 



d/2+k 

dz Cot d+1 (z) 

sphere Jd/2-k 



(_-]\d rd/2+k 
= K —^~ dzB d d+1 \ Z )7T COt (7TZ) 

a - Jd/2-k 

= f dz ^~ 1 J, d+1 {B d d+1) (d/2 + z)- B d d+1) (d/2 -z))tt taa(irz) 
Jo d\ 

= dz ^~^ 7T2; tan ttz JJ (z 2 - (j - 1/2) 2 ) 

f k 

= / dz P{d, z) tv tan7T2 , 
Jo 

(22) 

which is expressed in terms of elementary functions, (cf Das and Dunne, [17]. 
Kamela and Burgess, [18]). The even d result cannot be expressed in terms of 
multiple trig functions. (But see the section 9.) 



8. Negative modes and lattice considerations 



Expression (22) is not suitable for numerical computation. It is better to use 
a form involving the ip function. However, it is instructive when considering the 
behaviour as k increases, because for k > d/2 (allowed for odd d) the integrand 
diverges at z = d/2 + n, (n = 0, 1, . . .) and there appears to be a problem in the 
value of the integral. This is solved by taking the integral to be a contour one. 
Further, since the residues at the poles are integers, as the contour is changed, 
the left-hand side is uncertain only up to an integer multiple of 2iri so that the 
determinant is well defined. This is the same argument used by Gangolli, [19], in 
connection with the Selberg (^-function. 

Furthermore, splitting up the integral into principle part and residue contribu- 
tions, the latter are integer multiples of Tii and count the number of negative modes 
that are encountered as k increases. (I am thinking of k as a continuous variable.) 
The principle part corresponds to the ^function constructed from the moduli of 
the eigenvalues and the residues give the sign factor, (— l) t)A<0 . 



11 



Looking back to the eigenvalues, (2), the condition for a negative eigenvalue is 
the union of the N and D conditions, 

d 

J2m l <j + l-d/2, N 

i=i 

< j + l-(d+l)/2 / 
d (23) 

J2m t <j-d/2, D 

i=\ 

<j-(d+l)/2, 

with < j < k — 1, for a given k, (j G Z). 

The problem posed by (23) to find all m i > is the same as that underlying the 
notion of the Ehrhart polynomial, [20] , which is to find the number of non-negative 
integer lattice points, m { G Z d , in the integer dilation of the rational polytope, 

d 

rriiUJi < 1 , 

i=i 

where the uJ i are positive integers. In (23), we are dealing with problems in a 
standard integer simplex in ^-dimensions: all the U3 i are one. This is an easy and 
much studied example. Here, the integer dilations, t, are given by the right-hand 
sides of (23) and the answers are Ehrhart polynomials, in t. They are the binomial 
coefficients, [20]. This agrees with the residues in (22), as mentioned, if we 

note that, 

c^'+i) =( d V)- 



9. The N-D determinant ratio in the even, subcritical case 



For even d, as said, subtracting the N and D hemisphere log determinants 
removes the anomaly. This difference is easily found from (15) to be, 



Z' d (0,k) 



N-D 



= log 
= log 
= log 



r d (d/2 - k) 
r d (d/2 + k) 

T d+1 (d/2-k) T d+1 (d/2 + k + l) 
T d+1 {d/2 + k) F d+1 (d/2-k+l) 
Sm d+1 (d/2 + k) 



Sm d+1 (d/2-k) 



d even , 
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i.e. exactly the same form as the odd <i-sphere log determinant. 

The hemisphere approach puts a geometrical slant on Kurokawa's construction 
of the multiple sine function. 

The explicit calculation proceeds almost identically to the odd case, 



Z' d (0,k) 



d/2+k 

dz Cot d+1 (z) 

N-D Jd/2-k 

/ q \d r d/2+k 



d\ 



pa/ z^k 

/ dzBf +1 \z) 7T COt (TVZ) 

Jd/2-k 

= £ dz (B d d+1) (d/2 + z)- (d/2 -z))n cot (ttz) 

rk (-\)d d/2 ~ 1 

= / dz-r —7 nz cot nz 111 (z 2 -f), k < d/2 . 

Jo v-iy. f = \ y 

(24) 

Extending the range of k beyond its limit, d/2, by taking the integral as a 
contour one, still makes sense. This time, the pole residues count the difference 
in number of Neumann and Dirichlet negative modes which equals the number 
of Neumann negative modes in one dimension less, corresponding to the Ehrhart 
polynomial, (<£}«) = (?) -(*«-). 



10. Holographic aspect 

The integrand in the odd d expression, (22), is recognised as (proportional 
to) the continuation of the hyperbolic Plancherel measure on H d+1 , e.g. [21]. This 
equation is related to, but not the same as, a theorem of Kurokawa's, [22,23], in 
the special case of the trivial factoring of the (d + l)-dimensional hyperbolic bulk, 
H d+1 . The integral over the Plancherel measure came from the original functional 
relation of the Selberg ^-function which was then connected with the determinant 
of a 'linear' pseudo-differential operator constructed from the Laplacian on the 
dual space, S d+1 . (Sarnak gave the lowest dimension treatment in [24].) My result, 
on the other hand, has a holographic aspect as it relates a determinant on an S d 
boundary to quantities on the hyperbolic bulk. 

This can be made more explicit by recognising in (19) and (21) the ratio of two 

1/2 

determinants of the (d + 1) -dimensional spherical operators, {Yd+i + 1/4) ± k, 
computed again by considering the union of N and D (d + l)-hemisphere spectra, 
(m.aj + a D ± k) U (m.aj> + a N ± k), where a N = d/2 and a D = d/2 + 1. This quickly 
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gives, 

-logdetP 2fc = z;i(0,fc) 

d— sphere 

det[(y d+1 + l/4) 1/2 + A;l / (25) 

= lo § l ^ J±1 77172 7 + 2M ^ k ) ' 

det[{Y d+1 + l/4) 1/2 -k] 

as a holographic-type statement, the top line referring to <i-dimensional quantities 
and the bottom line to (d + l)-dimensional, apart from the polynomial, M(<i, k). If 
d is odd, there is no M term in (25). 

For odd d, the (d + l)-dimensional, purely bulk, formula that follows on com- 
bining (25) and (22), i.e. 

i det[(y d+1 + l/4) 1/2 + fc] f k 

log — — — jt -= / dzP(d, z) 7T tan nz (26) 

det[{Y d+1 + l/4) 1/2 -k] Jo 

is effectively given in [22], and written out in [25], Corollary 4.6. For Riemann 
surfaces, Cartier and Voros, [26], give equivalent equations. 

It is clear from (26), say by differentiating with respect to k to give a resolvent, 
that the residues at the poles of the integrand at n + d/2 are the degeneracies of 

1/2 

the eigenvalues, n + d/2, of (Y d+1 + 1/4) . Since it is a general theorem (see e.g. 
Camporesi, [21] §10, Helgason, [27]) that the residues of the Plancherel measure 
are (proportional to) the degeneracies on the dual compact space, here S d+1 , this 
proves that the polynomial P(d, z) is the Plancherel polynomial up to a factor, as it 
has here been computed to be, cf [23]. Also, the N U D hemisphere mode structure 
is equivalent to Theorem 3 in [23], for SO(<i + 2). 

Different bulks give different Plancherel measures (e.g. [23] and [28]). This is 
clear for the Yamabe case, k = 1, when the calculation for the other symmetric 
spaces can be followed through. Minakshisundaram, [29], gave an early analysis 
of the ^function for the unitary sphere and the other cases can be easily worked 
out. Cahn and Wolf [30], construct the heat-kernels and Ikeda [31] goes through 
the corresponding details of the (-functions. See also [23]. General and specific 
information can be found in the rather detailed Camporesi and Higuchi, [32], and 
[21]. 

For k 7^ 1, the analysis remains to be done. The equation equivalent to (26) 
obviously holds in terms of the Plancherel measure on the symmetric space and 
the Laplace-Beltrami operator on its compact dual. The question is deriving the 
identity (25) for a corresponding boundary GJMS operator, P 2k , given, possibly, as 
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a group representation intertwinor, cf Fontana, Branson and Morpurgo, [33]. One 
could conjecture a product structure similar to (1). 

The lattice considerations of section 8 again show the relation to the dimensions 
of the spherical representations of SO(d + 2). I also conjecture that this discussion 
extends to the other Plancherel measures and leads to Ehrhart polynomials for the 
various Cartan-Weyl weight polytopes, cf Welleda Beldoni, Beck and Cochet, [34]. 
Indeed the results of Panyushev, [35], Conway and Sloane, [36], and of Bachet, de 
la Harpe and Venkov, [37], can be used to obtain Theorem 4 of Kurokawa, [23], or, 
perhaps, vice versa. 

It is also possible to finitely factor the bulk. Apparently, some factorings cor- 
respond to spectral problems on the boundary that cannot be resolved classically, 
[6]. 



11. A generalisation 



It is clear that the formal developments do not depend on the particular values 
of the aij constants in (3) giving the GJMS operator. One can therefore define a 
more general (-function, Y d , from Z d , with any a ■, by, 

T d (a,k,a) 



Z' d (0,a,k,a) = log 



Qd(k) 



maintaining Barnesian normalisation with g v the necessary modulus. 

The N and D anomalies have opposite signs for odd d and the same signs 
for even d. To eliminate the anomaly, therefore, one can add the Z'ifS) for odd 
dimensions and subtract them for even, cf [12]. This motivates the combination, 

(-i) d 



'T d (a,k,a)~ 


-l 


"T d (a*, k,a)' 


. Qd(k) _ 







fc-l 

n 

3=0 



r d (a + a)r d (a-a) 



r d (a* + a,)r d (a*-a,) 



Pi 



(-i)° 



(27) 



which equality follows from (15). Here a and a* = d—a are dual boundary condition 
constants, such as a N and a D . 

The definition of the Barnes multiple sine function (see the appendix) suggests 
that the left-hand side of (27) be defined as the multiple sine function, &in d , for 
the Z d (-function, and then this equation is more compactly written, 

fc-i 

©in d (a, k, a) = J J Sin d (a + a^) Sin d (a — a^) . 
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12. Comments and conclusion 

The conformal anomaly and functional determinant of the GJMS operators, 
P 2 fc? ° n the sphere have been evaluated in explicit formal and analytical forms. As 
such, they should be useful as test cases for more general situations and also as 
fiducial expressions. 

Equation (22) agrees with the form extracted from the results of Guillarmou, 
[6], Theorems 1.3 and 4.2, derived on the basis that the _P 2 fc operators are the 
residues at a pole of the scattering operator on a hyperbolic bulk, of which the 
sphere is the conformal infinity. I have not needed to employ the Selberg ^-function. 

Equation (22) is also consonant with the expression derived by Diaz and Dorn, 
[7] in a boundary conformal field theory evaluation using dimensional regularisa- 
tion, with k now playing the role of conformal field dimension. These authors also 
produce an expression when d is even which should agree with the one here, up to 
a polynomial in k. 

Some history of ordinary sphere determinants was attempted in [38] and an 
improved technique, over that in [11], was given, but I have not used it. To the 
references in this paper, I should add those to Kurokawa mentioned above. 

The (^-function, (3), is an example of one defined by Friedman and Ruijsenhaas, 
[39], as a generalisation of Shintani's (^-function, [40]. The methods in the present 
work can be thought of as elementary and explicit derivations of some of the results 
given in this reference, applied to the GJMS operators. My approach is different 
and, in its domain of validity, simpler. 9 

Acknowledgments 
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Appendix 

In this appendix I place some, more or less standard, technical information 
on the Barnes functions that is best separated from the main text. I use Barnes' 
original definitions. 

9 The method of expanding in a can be applied to a general product of linear factors and will be 
described at another time. 
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The multiple ^functions are defined in terms Barnes' multiple T-function, 

^ P) (a) = ^logr». (28) 

I set ifj d (a) = ip^\a) for the first ^-function. Barnes, [13] §52 gives the recur- 
sion, 

*P d+1 (a + 1) = — ^ d (a) - - d CM a) (29) 



with 



\d 



( d (0,a)= { —^B d d \a). (30) 

I have put all parameters equal to 1 and rearranged the formula a little. 

Also, the recursion formula for the Barnes ^-function leads to the recursion for 
the multiple gamma function, [13], 

T d (a) _ T d+l (a) 



p d T d+1 (a + l)' 

which gives 

ip d+1 {a + 1) = Vd+i(a) - ^d(a) • 

Hence ^ 

ip d+1 (a) = — ^— ip d (a) - - C d (0, a) , (31) 

which can be iterated down to the standard function, ip. The result is, 

Ua) = («) ^(a) - £ ^ B^iC* - a) B<* >(a) , (32) 

where some products have been rewritten as Bernoulli polynomials. Pochhammer 
symbols could be used instead. 
The first three cases are, 

il> 2 (z) = -(z-l)1>(z) + z-- 

, , , z(z-l) , 18z 2 + 42z-17 
w) = g ^ 24 

, . , z(z-l)(z-2) (/ , 22z 3 -114z 2 + 167z -60 
M*) = ^ " ^(*) + n ' 

The general expression is easily programmed. 
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Similar expressions can be found in Kanemitsu et al, [41] and Onodera, [42]. 
Because of different normalisations and definitions, the formulae in some of these ref- 
erences are complicated by the appearance of derivatives of the Hurwitz ^-function. 

The formula (31) is useful for the evaluation of logr d (z) needed in the con- 
struction of the spherical log determinants. Thus 



r d (*i) 



(-i)d-i r*2 , n ( 33 ) 
/ dz(B^z)rl,(z) + Q d (z)) 



where the polynomial Q is given by 



= -(-l)^ 1 £ ^ B£&V - z) B<£\z) 
k=i 

and there are no poles of ip(z) in the integration range. 

Combinations of the r or ip functions often occur, in particular those leading 
to multiple trigonometric functions. These are usually defined using a differently 
normalised T-function. I have been using Barnes' original definition which is, 

log ^ = CM, 
Pd 

where p d is a normalising modulus. In [22], the T-function is defined without this 
factor, 

Pd 

which is Vardi's notation, [43]. 10 In terms of G d , the multiple sine function is 
defined by, 

Si n dW ^-» < - 1,d 



GJz) 

and I note that for even d this is independent of p d . 

The ^-function, (28) , also does not depend on p d and one can define a multiple 
cotangent function by, 

Cot d (z) = ^logSm d (z) {u) 
= -*l> d {z)-{-l) d *l> d {d-z), 



10 This topic is beset by conflicting definitions and notations. The G d should not be confused 
with the quantity having the same symbol in Vigneras, [44]. Vardi denotes this quantity by T d . 
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which has the nice property that it satisfies the homogeneous recursion, 

Cot d+1 (a) = ^^Cot d (a), 
easily proved from the ip recursions using the symmetry, 

Bf\d-a) = {-l) d Bf\a). 
Hence one obtains the explicit connection with the ordinary cotangent, 

C ° td(a) = ^^BdlM^cot^a). (35) 

These and other relations can be found in [25], and elsewhere, derived in various 
ways. For example, from the definition (34) one has the duality (or complementar- 
ity) symmetry, 

Cot d (d - z) = (-l) d Cot d (z) . 

References. 

1. Graham,C.R., Jenne,R., Mason,L.J. and Sparling, G.A.J. J. Lond. Math. Soc. 

46 (1002) 557. 

2. Juhl,A. Families of Conformally Covariant Differential Operators, Q-Curvature 

and Holography, (Birkhauser, Basel, 2009). 

3. Gover,A.R. Laplacian operators and Q-curvature on conformally Einstein man- 

ifolds ArXiv:math.DG/0506037. 

4. Diaz,D.E. JHEP 7 (2008)103. 

5. Branson,T.P. Trans. Am. Math. Soc. 347 (1995) 3671. 

6. Guillarmou,C. Am. J. Math. 131 (2009) 1359. 

7. Diaz,D.E. and Dorn,H. JHEP 0705 (2007) 46. 

8. Graham,C.R. SIGMA 3 (2007) 121. 

9. Morpurgo,C. Duke Math. J. 114 (2002) 477. 

10. Chang,P. and Dowker,J.S. Nucl. Phys. B395 (1993) 407. 

11. Dowker,J.S. Comm. Math. Phys. 162 (1994) 633. 

12. Dowker,J.S. J. Math. Phys. 35 (1994) 4989. 

13. Barnes,E.W. Trans. Camb. Phil. Soc. 19 (1903) 374. 

14. Dowker.J.S. Entanglement entropy for even spheres. ArXiv: 1009.3854. 

15. Norlund,N.E. Acta Mathematica 43 (1922) 121. 



19 



16. Barnes,E.W. Trans. Camb. Phil. Soc. 19 (1903) 426. 

17. Das,A. and Dunne,G.V. Phys. Rev. D74 (2006) 044029. 

18. Kamela,M. and Burgess,C.P. Can. J. Phys. 77 (1999) 85. 

19. Gangolli,R.///mozs J. Math. 21 (1977) 403. 

20. Beck,M. and Robins, S. Computing the Continuous Discretely, (Springer, New 

York, 2007). 

21. Camporesi,R. Phys. Rep. 196 (1990) 1. 

22. Kurokawa,N. Proc. Jap. Acad. A67 (1991) 61. 

23. Kurokawa,N. Proc. Jap. Acad. A68 (1992) 256. 

24. Sarnak,P. Comm. Math. Phys. 110 (1987) 113. 

25. Kurokawa,N. and Koyama,S. Forum. Math. 15 (2003) 839. 

26. Cartier,P. and Voros,A. Grothendieck Festschrift, v. 2, p.l-69,(Birkhauser,1991). 

27. Helgason,S. Asterisque hors serie, 1985, 151. 

28. Wakayama,M. Hiroshima Math. J. 15 (1985) 260. 

29. Minakshisundaram.S. Can. J. Math. 4 (1952) 26. 

30. Cahn,R.S. and Wolf, J. A. Comm. Mat. Helv. 51 (1976) 1. 

31. Ikeda,A. Kodai Math. J. 18 (1995) 57. 

32. Camporesi,R. and Higuchi,A. J.Geom. and Physics 15 (1994) 57. 

33. Fontana,L., Branson,T.P. and Morpurgo,C. Moser-Trudinger and Beckner- 

Onofri inequalities on the CR sphere, ArXiv:math.AP 0712/3905. 

34. Welleda Beldoni,M., Beck,M. and Cochet, C. Discrete and Comp.Geom. 35 

(2006) 55. 

35. Panyuchev,P.I. Func.Anal.Appl. 28(1994)293. 

36. Conway,J.H. and Sloane,N.J.A. Proc. Roy. Soc. A453 (1997) 2369. 

37. Bacher,R., de la Harpe,P. and Venkov,B. Ann. I'inst. Fourier 49(1999)727. 

38. Dowker,J.S. and Kirsten, K. Analysis and Appl. 3 (2005) 45. 

39. Friedman,E. and Ruijsenhaars,S. Adv. in Math. 187 (2004) 362. 

40. Shintani,T. J. Fac.Sci.Univ. Tokyo 24 (1977) 167. 

41. Kanemitsu,S., Kumagai,H. and Yoshimoto,M. The Ramanujan J. 5(2001)5. 

42. Onodera,K. Adv. in Math. 224 (2010) 895. 

43. Vardi,I. SI AM J. Math. And. 19 (1988) 493. 

44. Vigneras,M.F. Asterisque 61 (1979) 235. 



20 



